function prof = CAT_prof1(b,X,m,c)

prof = mean(m*(X*b-c).*((X*b > c)-1));